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It  is  projected  that,  within  the  next  few  years,  China  will 
launch  its  first  natural  resource  satellite  (ZY-1  satellite) . 
Satellite  and  ground  reception  as  well  as  data  processing 
application  systems  also  follow  it,  beginning  preparations  for 
construction.  Satellite  data  preprocessing  is  an  important 
subsystem  in  ground  application  systems.  Preprocessing  technology 
is  basically  weak  in  China.  Beginning  the  development  of  research 
associated  with  preprocessing  geometrical  correction  algorithms  for 
ZY-1  satellite  products  is  very  necessary. 

I.  ZY-1  SATELLITE  IMAGERY  GEOMETRICAL  CHARACTERISTICS 

ZY-1  satellites  operate  in  near  circular  solar  synchronc  is 
orbits.  Their  main  useful  load  is  CCD  cameras  and  IR-MSS  scanning 
devices.  They  are  capable  of  separately  obtaining  multiple 
spectrum  imagery  of  the  earth.  CCD  cameras  also  possess 
panchromatic  wave  band  imagery.  Its  various  wave  bands  are  all 
6000  image  point.  In  vertical  observations,  ground  resolutions  are 
19.4m.  IR-MSS  scanning  devices  are  1536  image  point.  Ground 
resolution  is  78m.  CCD  camera  directional  lenses  permit  side  views 
on  the  two  sides  of  orbits  deviating  from  the  nadir.  In  space 
perpendicular  to  orbital  direction — using  2 '  interval  partially 
blocked  slopes,  maximum  side  view  angles  are  ±32°.  ZY-1  satellite 
control  system  attitude  assurance  error  <0.15°. 

CCD  camera  form  is  central  projection.  At  a  certain  instant 
(that  is,  x=0)  the  expansion  derived  for  form  equations  is: 
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In  this,  f — camera  primary  length;  x,y — image  coordinates  in 
sensor  coordinate  system;  all,  al2,  . ..,  a33 — attitude  angle 

rotation 

matrix  coefficients;  — ground  coordinates  of 

target  points; 

xry,j,  — coordinates  of  points  under  satellite. 


IR-MSS  infrared  scanning  instrument  geometrical  form  is  a 
panoramic  projection  as  time  x=0.  The  corresponding  structural 
equations  are: 
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In  this,  0 — scanning  angle. 


II.  ZY-1  SATELLITE  IMAGERY  DEFORMATION  ERRORS 

ZY-1  satellite  imagery  deformations  can  be  roughly  divided 
into  two  classes — internal  errors  and  external  errors.  Interior 
errors  are  given  rise  to  by  the  structural  characteristics  of/53 
sensor  devices  becoming  nonideal  or  their  indices  deviating  from 
nomimal  numerical  values — for  example,  standardized  primary  camera 
lengths  being  different  from  the  actual  ones;  noncontinuity  of 
surface  features  on  imagery  given  rise  to  by  bidirectional  scanning 
errors  associated  with  multiple  spectrum  scanning  devices;  the 
existence  of  such  phenomena  as  intersections,  omissions,  as  well  as 
overlaps  between  image  elements  and  scanning  lines,  and  so  on,  and 
so  on.  Among  these  errors,  systemic  ones  should  go  through  sensor 
device  inspections  in  order  to  be  calibrated  or  be  given  correction 
formulae  to  supply  to  processing  systems  for  corrections.  As  far 
as  external  errors  are  concerned,  that  is,  sensor  device  errors 
with  respect  to  positioning  during  imagery  processes,  they 
primarily  include  the  several  types  below. 

1.  Changes  in  Sensor  Device  External  Azimuth  Elements 

With  regard  to  the  analysis  of  the  deformation  errors  in 
question,  it  is  possible  to  carry  out  differentiation  in  respect  to 
the  two  ends  of  collinear  sensor  equations  by  directional 
displacements  associated  with  the  coordinates  x,y,z  of  points  under 
satellites  and  the  three  sensor  attitude  angles  i|»  ,  .w,  k  being 
variables.  One  obtains: 
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Again,  on  the  basis  of  sensor  structural  equation  rotation 
matrices,  it  is  then  possible  to  derive  the  image  point 
displacement  formulae  given  rise  to  by  dx,  dy,  dz,  dt|r,  dia,  and  dK. 

2 .  Imagary  Errors  Given  Rise  to  by  the  Spin  of  the  Earth 
When  ZY-1  satellites  make  scanning  movements  from  north  to 
south,  due  to  the  earth,  there  is  also  a  spinning  from  west  to  east 
at  the  same  time.  As  a  result,  the  loci  of  points  under  the 
satllites  gives  rise  to  a  bending  and  sloping  toward  the  west. 
Simultaneously,  because  the  geographical  latitudes  associated,  with 
points  under  satellites  are  different,  the  linear  velocities 
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associated  with  the  spinning  earth  VE  at  the  locations  in  question 
are  also  different.  As  a  result,  errors  given  rise  to  are  capable 
of  going  through  corrections,  solving  for  transverse  displacement 
values  dye  of  any  scanning  line  in  an  image  relative  to  the  first 
scanning  line  and  the  average  image  deflection  angle  0y  . 


3 .  Image  Deformations  Given  Rise  to  by  the  Curvature  of  the 
Earth 

With  regard  to  the  influence  of  the  curvature  of  the  earth,  it 
is  possible  to  make  use  of  height  differences  Ah  given  rise  to  by 
the  earth's  degree  of  curvature,  which  exist  for  images  projected 
on  a  spherical  reference  system,  or  to  make  use  of  points  on  the 
ground  relative  to  projection  planes,  deriving  imagery  earth 
curvature  image  point  displacement  formulae  and  applying 
corrections. 

III.  ZY-1  SATELLITE  IMAGERY  GEOMETRICAL  CORRECTIONS 

As  far  as  ZY-1  satellites  opting  for  the  use  of  digital 
imagery  is  concerned — with  regard  to  imagery  geometry  calibrations 
carried  out  in  association  with  imagery  errors  given  rise  to 
because  of  satellite  orbital  deviations  and  sensor  attitude 
changes--on  the  basis  of  digital  correction  models,  it  is  possible 
to  distinguish  "reference  methods"  and  "nonreference  methods". 
Parameter  methods  primarily  refer  to  col inear  equation  methods. 
Nonreference  methods  mostly  indicate  multiple  term  type  correction 
methods  as  well  as  random  field  interpolation  value  methods,  and  so 
on.  General  processes  associated  with  digital  correction  methods 
are  as  shown  in  Fig.l. 
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Fig.l  (l)  Input  Original  Pattern  (2)  Set  Up  Correction  Model 
(3)  Precisely  Determine  Output  Pattern  Range  (5)  Correction  of 
Image  (Element  Geometrical  Positions  (6)  Resampling  of  Pattern 
Greynejss  Values 

Satellite  imagery  geometrical  correction  preprocessing 
normally  is  also  divided  into  rough  processing  and  precision 
processing.  Rough  processing  generally  indicates  making  use  of  a 
number  of  geometrical  parameters  which  can  be  predicted  during 
imagery  formation  processes,  associated  with  scanning  images,  in 
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order  to  directly  construct  correction  formulae  and  carry  out 
calibrations.  As  far  as  comparisions  between  precision  processing 
and  rough  processing  are  concerned,  the  differences  are  that  /54 
precision  processing  requires  making  use  of  ground  control  points 
to  act  as  correction  bases.  On  the  basis  of  requirements 
associated  with  the  "Overall  Plan  for  Natural  Resources  No.l 
Satellite  Application  System  Technology",  ZY-1  satellite  imagery 
products  are  divided  into  6  levels: 


0  Level: 

1  Level: 

2  Level r.  - 

3  Level: 

4  Level: 

5  Level: 


Divided  amplitude  data  products  which  have  not 
undergone  any  calibration 

Products  Having  Undergone  Radiation  Calibration 
System  Level  Geometrical  Calibration  Products 
Two  Dimensional  Geometrical  Calibration  Products 
Using  a  Series  of  Ground  Control  Points  as  a 
Foundation 

Products  Making  Use  of  Ground  Elevation  Models  to 
Act  as  Geometrical  Corrections 
High  Level  Nonstandard  Products 


The  large  amounts  of  ZY-1  satellite  preprocessing  products 
should  be  divided  into  the  two  classes  of  imagery  products. 
Besides  this,  based  on  user  requirements,  use  is  made  of  small 
numbers  of  ground  control  points  to  act  as  a  third  class  of  batch 
precision  processing  products.  With  regard  to  more  advanced 
precision  processing  products,  users  need  not  directly  make  use  of 
such  supplementary  data  as  satellite  orbits,  attitudes,  and  so  on, 
using  polynomial  equations  for  autonomous  processing  or  requiring 
special  processing  by  processing  centers. 

Level  2  ZY-1  satellite  products  directly  make  use  of  satellite 
orbits,  sensor  attitude  data,  and  structural  equations  based  on 
sensors  to  construct  geometrical  correction  formulae  to  carry  out 
corrections.  As  far  as  the  actual  situations  given  rise  to  by  the 
combined  factors  of  unknown  attitude  errors  and  errors  associated 
with  orbits  announced  in  advance  are  concerned,  mean  square 
deviations  associated  with  predeterimed  views  should  be  within  a 
range  of  800m. 

With  respect  to  ZY-1  satellite  imagery  level  2  products,  the 
methods  of  geometrical  correction  based  on  different  sensors  are  as 
follows. 

1.  IR-MSS  Scanning  Device  Imagery  Product  Rough  Processing 

Based  on  experience  associated  with  the  U.S.  Landset 
satellite,  ZY-1  satellite  imagery  is  capable  of  separately 
supplying  users — after  opting  for  the  use  of  two  types  of 
correction  methods — with  level  2  products. 

(1)  Slope  Elimination  Corrected  Products 

To  primarily  carry  out  corrections  of  errors  associated  with- 
imagery  distortion,  single  event  expansions  and  contractions,  and 
deflection  deformations  given  rise  to  by  such  factors  as  the  spin 
of  the  earth,  image  element  x  and  y  direction  dimensions  not  being 
equal,  as  well  as  satellite  orbital  plane  deflection  angles  (that 
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is,  the  angles  included  between  the  axses  of  spin  of  the  earth) , 
and  so  on,  the  formula  is: 


0  cosa  Sina 

M  —  Sina  cosa 


In  this,  x,y — original  imagery  coordinates  associated  with  a 
certain  image  element;  x',y' — corrected  imagery  coordinates 
associated  with  the  corresponding  image  element;  My — ZY-1 
satellite  imagery  aspect  ratio;  a — satellite  movement  direction 
angle;  0y — mean  imagery  distortion  angle;  Mo — rotation  matrix 
associated  with  imagery  top  level  after  making  corrections. 

(2)  Map  Projection  Correction  Methods 

The  principle  is  that  one  first  defines  an  image  grid  on  the 
original  imagery.  Making  use  of  sensor  device  movement  parameters 
that  can  be  predicted,  one  uses  grid  point  coordinates,  converting 
them  into  geographical  latitude  and  longitude  in  accordance  with 
colinear  equations.  Going  a  step  further,  they  are  converted  to 
map  projection  system  coordinates,  using  grid  lattice  points  to  act 
as  artificial  control  points  in  order  to  calculate  a  set  of 
corrected  polynomial  coefficients.  Then,  calibrations  are  carried 
out  on  the  basis  of  polynomials  for  all  points  with  respect  to  the 
imagery.  As  far  as  concrete  algorithms  are  concerned,  it  is 
possible  to  opt  for  the  use  of  simplified,  high  speed  algorithms. 

2 .  Rough  Processing  of  CCD  Camera  Level  2  Product 

Defining  satellite  height  changes  as  dzs,  define  photography 
station  translations  given  rise  to  by  satellite  speed  changes  as 
dxs  and  dys»  Defining  camera  directional  lens  positioning  errors 
as  dty  and  dio,  it  is  possible  to  construct  image  point  displacement 
formulae,  carrying  out  corrections  with  respect  to  errors  given 
rise  to  by  external  azimuth  elements. 

With  regard  to  the  influences  of  the  curvature  of  the  earth, 
it  is  possible  to  take  imagery  and  project  it  onto  spherical  shaped 
reference  systems  and  achieve  corrections. 

The  influences  of  the  spin  of  the  earth  primarily  create 
imagery  scanning  line  transverse  displacements  in  the  y  direction 
and  scale  enlargements  in  the  x  direction.  As  far  as  translation 
improvements  associated  with  ZY-1  satellite  imagery  are  /55 
concerned,  it  is  possible  to  opt  for  the  use  of  line  by  line 
scanning  improvement  methods  in  the  same  way  as  SPOT  satellite 
imagery . 

In  regard  to  the  contents  of  ZY-1  satellite  CCD  imagery  level 
2  products,  the  effects  of  such  things  as  the  undulations  of 
terrain,  sensor  device  attitude  changes,  angles  of  drift,  and  so 
on,  are  preserved.  Although  geometrical  precision  is  relatively 
low,  terrain  change  data  is  retained,  however.  This  is 
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advantageous  to  three  dimensional  observations  and  the  utilization 
of  original  terrain  data. 

With  respect  to  ZY-1  satellite  orbital  characteristics  and 
sensor  device  imagery  geometrical  configurations,  they  are 
respectively  similar  to  the  U.S.  Landset  satellite  and  the  French 
SPOT  satellite.  Therefore,  in  making  use  of  satellite  prediction 
parameters  to  carry  out  rough  processing  associated  with 
geometrical  preprocessing — in  terms  of  the  principles  of 
processing — they  are  the  same  as  the  foreign  satellites.  In 
computer  software  algorithms  and  software  writing,  design  will  be 
autonomous . 

As  far  as  geometrical  calibration  of  ZY-1  satellite  level  3 
product  is  concerned,  it  is  possible  to  make  use  of  simplified 
correction  models — for  example,  opting  for  the  use  of  planar  grid 
distortion  model  interpolation  methods  or  other  high  speed 
algorithms,  making  use  of  as  few  ground  control  points  or  ground 
control  point  imagery  banks  as  possible  to  carry  out  corrections. 


Fig. 2  (1)  Level  1  Product  or  Original  Data  (2)  Colinear 

Equation  Model  (3)  Precision  Processing  Level  3  Product  (4) 
Ground  Control  Point  (5)  Greyness  Level  Signature  Sampling 

With  respect  to  satellite  imagery  geometrical  corrections — 
after  conversion  of  geometrical  positions — there  is  generally  a 
need  in  all  cases  to  carry  out  sampling  of  imagery  greyness  values 
again.  During  previous  satellite  imagery  preprocessing — on  a 
foundation  of  roughly  calibrated  products — further  precision 
processing  is  generally  always  done.  For  example,  each  level  of 
product  associated  with  the  French  SPOT  satellite  imagery  is  a 
reprocessing  into  higher  level  products.  In  this  way,  precision 
processed  products  will  all  go  through  several  interations  of 
greyness  value  resampling.  Although  option  is  made  for  the  use  of 
interpolation  prinicples  in  order  to  reduce  greyness  level  drops, 
due  to  multiple  iterations  of  greyness  resampling,  however,  imagery- 
greyness  levels  are  still  created  which  must  necessarily  go  down. 

In  order  to  escape  from  this,  the  three  levels  of  ZY-1 
satellite  products  should — from  satellite  level  1  imagery  products 
and  not  making  use  of  level  2  products — go  through  calibration 
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models  and  ground  control  points.  In  conjunction  with  this, 
methods  used  produce  precision  processing  products  one  time, 
method  flow  chart  is  seen  in  Fig. 2. 


The 
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